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Abstract. The atmospheric parameters and iron abundance of the Sloan Digital Sky Survey (SDSS) spectrophotometric standard 
star BD +17 4708 are critically examined using up-to-date Kurucz model atmospheres, LTE line formation calculations, and 
reliable atomic data. We find r eff = 6141 ± 50 K, logg = 3.87 ± 0.08, and [Fe/H] = -1.74 ± 0.09. The line-of-sight interstellar 
reddening, bolometric flux, limb-darkened angular diameter, stellar mass, and the abundances of Mg, Si, and Ca are also 
obtained: E(B - V) = 0.010 ± 0.003, f bol = (4.89 ± 0.10) x 10"' erg cirT 2 s~\ 6 = 0.1016 ± 0.0023 mas, M = 0.91 ± O.O6M , 
[Mg/Fe] = 0.40 ± 0.10, [Si/Fe] = 0.35 ±0.11, [Ca/Fe] = 0.36 ±0.11. This star is a unique example of a moderately metal- 
poor star for which the effective temperature {T^g ) can be accurately constrained from the observed spectral energy distribution 
(corrected for reddening). Such analysis leads to a value that is higher than most spectroscopic results previously reported in the 
literature (~5950 K). Interstellar reddening was estimated using various prescriptions, including an analysis of interstellar lines. 
The surface gravity of the star was inferred from the fitting of the wings of the Mgib lines. We used transition probabilities 
measured in the laboratory and reliable damping constants for unblended Fe lines to derive the iron abundance using both Fe i 
and Fen lines. We find that the ionization balance of Fe lines is satisfied only if a low r eff (~ 5950 K) is adopted. The mean 
iron abundance we obtain from the Fen lines corresponds to Ap c = 5.77 ± 0.09 ([Fe/H] = -1.74 for our derived Af c ,o = 7.51) 
while that from the Fei lines is Ap e = 5.92 ± 0.11, and therefore with our preferred T c h (6141 K), the discrepancy between 
the mean iron abundance from Fe i and Fe n lines cannot be explained by overionization by UV photons as the main non-LTE 
effect. Interestingly, the Fe i excitation balance is satisfied with a r c ff only slightly warmer than our preferred solution and not 
with the lower value of 5950 K. We also comment on non-LTE effects and the importance of inelastic collisions with neutral H 
atoms in the determination of oxygen abundances in metal-poor stars from the 7774 A O i triplet. 
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1. Introduction 

The derivation of stellar chemical compositions relies on the 
accurate determination of the atmospheric parameters T e g (ef- 
fective temperature) and logg (surface gravity). These quanti- 
ties may be inferred either from the stellar spectrum or by semi- 
empirical methods that are normally based on photometric and 
astro me trie measurements. Often, the photometric and astro- 
metric parameters are used as first approximations and the final 
solution is found iteratively with the help of the spectrum. This 
tuning of parameters is, however, model-dependent, and may 
lead to erroneous conclusions if the models are inadequate. 

Most abundance analyses of FGK stars are made using ho- 
mogeneous plane-parallel model atmospheres and LTE (local 
thermodynamic equilibrium) line formation. However, recent 
abundance analyses using hydrodynamical model atmospheres 
and non-LTE line formation have demonstrated that the effects 
of surface inhomogeneities and departures from LTE on abun- 
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dance analyses are not negligible in the Sun and solar-type stars 
of different metallicities (e.g., Asplund & Garcia Perez 2001, 
Korn et al. 2003, Allende Prieto et al. 2004b, Asplund 2005), 
being probably dramatic for very metal-poor stars (Shchukina 
etal.2005). 

The F8-type star BD +17 4708 has been chosen as a spec- 
trophotometric standard, either primary or secondary, for vari- 
ous systems (e.g., Oke & Gunn 1983, Rufener & Nicolet 1988, 
J0rgensen 1994, Zhou et al. 2001). In particular, this star is the 
primary standard of the Sloan Digital Sky Survey (SDSS) pho- 
tometric system (Fukugita et al. 1996; Gunn et al. 1998, 2006; 
Smith et al. 2002). BD +17 4708 is one of the few stars, and 
the only subdwarf, with very accurate absolute fluxes (Bohlin 
& Gilliland 2004b), which allows us to test, in an independent 
way, the models and different atmospheric parameters that have 
been derived for it. In particular, it is a unique case in which the 
effective temperature may be well constrained. Also, it can be 
quite useful to know, with high accuracy, what the fundamen- 
tal parameters and overall chemical composition of this star 
are, given that its model atmosphere and predicted fluxes may 
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Table 1. Literature data for BD +17 4708. Note that there are 
two sets of parameters given by Boesgaard et al. (1999): those 
given in the King (1993, K93) scale, and those given in the 
Carney (1983a,b; C83) scale. In this Table, only values derived 
for the first time (although mostly taken from the more recent 
papers) by each author(s) are shown to avoid duplicity. 



be used to complement the observed spectral energy distribu- 
tion in the transformation of observed magnitudes into physical 
fluxes in the SDSS. 

BD +17 4708 has been studied by several groups (Table^, 
who have derived effective temperatures between 5800 K and 
6200 K, [Fe/H] values between -2.0 and -1.4 and a surface 
gravity of logg ^ 4.0 (Fig. Q. 1 In this paper, we present a 
detailed determination of the atmospheric parameters and the 
iron abundance of BD +17 4708. A critical comparison with 
values previously reported in the literature is provided in the 
Appendix. Due to their relevance for studies of stellar interiors, 
we also derive the abundances of Mg, Si, and Ca, and discuss 
the uncertainties in the determination of the oxygen abundance 
from the 7774 A triplet, in particular the importance of inelastic 
collisions with neutral H atoms in the non-LTE computations 
involved. Our analysis is based on the observed spectral energy 
distribution and a high resolution, high signal-to-noise (S/N) 
spectrum of the star. 

2. Fitting of the spectral energy distribution 

Bohlin & Gilliland (2004b, hereafter BL04b) have measured 
the spectral energy distribution of BD +17 4708, from the UV 
(0.17 yum) to the near IR (1.0 fim), with respect to the three pri- 
mary standards of the Space Telescope Imaging Spectrograph 
(STIS) on the Hubble Space Telescope (HST). The accuracy of 
BL04b measurements is better than 0.5% and therefore, in the 
absolute scale, their fluxes for BD +17 4708 are as accurate 
as those of the three white dwarfs used as primary standards. 
The fluxes measured with STIS for these three white dwarfs 



1 Throughout this paper, we use the standard definitions = 
log{N x IN y ) - log(N x /N Y )e and A x = \og(N x /N H ) + 12, where N x is 
the number density of the element X. The surface gravity g in log g is 
given in cgs units. 
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Fig. 1. Distribution of the r e fj, [Fe/H], and logg values found 
in the literature for the star BD +17 4708, as given in Tabled 



have uncertainties that range from 2.5% in the UV to 1% at 
longer wavelengths, according to BL04b. Thus, the spectral en- 
ergy distribution of BD +17 4708 is accurate at the level of 2% 
or even better in some spectral regions. 

In Fig. |2 the observed fluxes from BL04b are shown along 
with theoretical flux distributions computed by R. L. Kurucz 2 
after applying an E(B — V) = 0.01 reddening (see s32. II ac- 
cording to the Fitzpatrick (1999) parameterization with R v - 
A V /E(B - V) = 3.1. The models have logg = 3.87 and 
[Fe/H] = -1.74, values that we derive in ^3] The theoretical 
fluxes have been empirically scaled to the observed one using 
the reddest 1000 A, i.e., they were divided by the mean ratio 
of theoretical to observed fluxes from 9000 A to 10000 A in 
each case. The scaling factor, s, is directly related to the stellar 
angular diameter, 8, by s = 6 2 /4. The spectrum from BL04b 
was smoothed to approximately match the resolving power of 
the Kurucz model fluxes. 

As shown in Fig. |2] up-to-date Kurucz models, which in- 
clude newly computed opacity distribution functions, an a- 
element enhancement consistent with the mean [a/Fe] ratio we 
obtain ( §I3.7> . and without convective overshooting, accurately 
reproduce (within 1%) the observed flux distribution at wave- 
lengths longer than 5000 A, as long as the correct T^ - 6141 K 
is adopted (^3- In the range 4000 A< A < 5000 A the 
model underestimates the flux by about 1.5% while at shorter 
wavelengths the fit is reasonable on average although several 
strong lines are not well fitted. Interestingly, a comparison of 
a MARCS model 3 with a Kurucz model of parameters near 
those of BD +17 4708 (the closest point in the grid we found 



2 Up-to-date models are available at http://kurucz.harvard.edu The 
characteristics of the models are explained in Kurucz (1970, 1979). 

3 Downloaded from http://marcs.astro.uu.se 



Ramirez et al.: The SDSS standard BD +17 4708 



3 



7.5 



6.0 



4.5 



3.0 



1 .5 



a 5.0 



3 
J* 



< 



0.0 



-5.0 



BD +17 4708 

[Fe/H]= 1.74 

log g = 3.87 
E(B-V)=0.0i 
T.„ (K) 

6141 [no] 

5950 [noj 

6222 [ov] 






i — i i i i | i i i i | i — i — | — i — i — i — i — | — t — i — i — i — = 


:, i , , , , i , >W 


^7^+^^ 



2000 3000 4000 5000 6000 7000 8000 9000 

MA) 

Fig. 2. Upper panel: The spectral energy distribution measured by Bohlin & Gilliland (2004b) is shown with the filled circles. 
Alpha-element enhanced ([a/Fe]=+0.4) no-overshoot [no] Kurucz's models of [Fe/H] = -1.74, logg = 3.87, reddened by 
E(B-V) = 0.01 according to the Fitzpatrick (1999) parameterization, and two effective temperatures: T e g = 6141 Kand5950K, 
are shown with the solid and dot-dashed lines, respectively. An overshoot [ov] Kurucz model of T e s = 6222 K is shown with the 
dotted line (log g and [Fe/H] for this model are slightly different to make them consistent with the higher T^ff). Lower panel: As 
in the upper panel for the difference (percentile) between observed and theoretical fluxes. The hydrogen lines have been excluded. 



is: r eff = 6000 K, logg = 4.0, [Fe/H] = -2.0) shows that 
the MARCS model predicts fluxes in this region that are larger 
by about 1%, which would reduce the difference somewhat. 
Kurucz and MARCS models seem to predict roughly the same 
fluxes everywhere else. 

The Kurucz overshoot model shown in Fig. [2] has the T e g 
that best fits the data for overshoot models. Clearly, it does 
not reproduce very well the observational data. The overshoot 
model does not have a-element enhancement but this has a 
much smaller effect on the shape of the spectral energy distri- 
bution compared to the switch to a no-overshoot model. Note 
also that the overshoot model that fits best the observations 
(that shown in Fig. O is hotter by about 80 K compared to the 
no-overshoot model, thus introducing a systematic error in the 
T e ff derived from these fits. For that reason, we prefer to adopt 
Kurucz no-overshoot models hereafter. Note also that adopting 
a Kurucz model of T e g = 5950 K, as suggested by previous 
spectroscopic studies (Tableland Fig.[0, results in a severe 
discrepancy with the observational data. 

Although it is not shown in Fig. [2 there is a slight degener- 
acy between T e s and E(B - V) in the model fits to the observed 
flux distribution. Roughly speaking, increasing r e ff is equiva- 
lent to decreasing E(B - V), as they both result in higher fluxes 
in the UV-blue regions while leaving the red and infrared fluxes 
nearly unchanged (after scaling). For instance, equally good fits 
to the data can be obtained with r e fr - 6050 K if no redden- 
ing is assumed or with r eff - 6150 K if E(B - V) = 0.010 is 



adopted. It is, thus, important to constrain the E(B - V) value 
independently. 

2.1. Reddening 

Since BD +17 4708 is at a distance of about 120 pc (its 
Hipparcos parallax is 8.43+1.42 mas), the E(B-V) value is ex- 
pected to be negligible or small. The Local Bubble, a region de- 
void of dense gas extends approximately 60 pc in the direction 
of BD +17 4708 (Lallement et al. 2003). At the Local Bubble 
boundary, a significant increase in dense interstellar medium 
(ISM) material is observed toward several stars in the general 
direction of BD +17 4708. This ISM material can be observed 
as narrow absorption lines in high resolution spectra, and de- 
pending on the strength of the absorption, it may be expected 
to cause a small, but measurable, amount of reddening. 

Interstellar gas in the line of sight of BD +17 4708 is 
evidenced by the interstellar Nai Dl (~ 5896 A) and D2 
(~ 5890 A) lines shown in Fig. Due to the high radial ve- 
locity of the star (-291 km s _1 ), the ISM absorption lines at 
-13 km s _1 are significantly displaced from the strong stellar 
features. Other than the two stellar lines and two interstellar 
lines, all of the remaining features seen in Fig.|3]are caused by 
telluric water vapor. 

Although the telluric H2O lines are relatively weak, they 
need to be modeled and removed from the spectrum, in order 
to obtain a high precision measurement of the Na 1 ISM column 
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density. A relatively simple model of terrestrial atmospheric 
transmission (AT - Atmospheric Transmission program, from 
Airhead Software, Boulder, CO) developed by Erich Grossman 
is used to fit the telluric water vapor lines. This forward model- 
ing technique to remove telluric line contamination in the vicin- 
ity of the Na i D lines is described in detail by Lallement et al. 
(1993), in which a more sophisticated terrestrial atmospheric 
model was employed. As can be seen in Fig. [3] the AT program 
is very successful at modeling the terrestrial absorption in the 
spectrum of BD +17 4708. Note that the stellar lines happen to 
fall in an area free of contaminating lines, and the two interstel- 
lar lines are only slightly blended with water vapor absorption 
in the wings of the ISM absorption. 

The interstellar lines found in the BD +17 4708 spec- 
trum were modeled using standard methods (see, e.g., §2.2 in 
Redfield & Linsky 2004a). A single Gaussian absorption com- 
ponent is fit to both Na i D lines simultaneously using atomic 
data from Morton (2003), and then convolved with the instru- 
mental line spread function. Fitting the lines simultaneously 
reduces the influence of systematic errors, such as continuum 
placement and contamination by weak features. The free pa- 
rameters are the central velocity (v), the line width or Doppler 
parameter (b), and most importantly, because it can be used 
to estimate the reddening along the line of sight, the column 
density (AO of Nai ions toward BD +17 4708. The best fit is 
shown in Fig. [3] where v = -13.315 + 0.031 km s _1 , b = 
3.922±0.052kms- 1 ,andlogJV(Nai) = 1 1.4776±0.0034cirT 2 . 
Due to the high S/N, systematic errors probably dominate over 
the statistical errors given above. 

The measured line width, a consequence of thermal and 
non-thermal, or microturbulent, broadening (see Redfield & 
Linsky 2004b), is wider than typically found for cold ISM 
clouds (Welty et al. 1994). This is likely due to unresolved 
interstellar components along the line of sight. Although this 
observation is considered to be high spectral resolution from a 
stellar perspective (A/6A ^ 60, 000, see 33. 1> . from an inter- 
stellar perspective it is moderate resolution because the narrow 
and closely spaced component structure typical of the ISM is 
best observed at much higher resolution (A/ 6 A ^ 500,000 - 
1 , 000, 000). Therefore, we are likely seeing blending of several 
ISM components along the line of sight toward BD +17 4708. 
In fact, 7r Aqr, a star in the same part of the sky as BD + 1 7 4708 
(A6 ~ 15°), but more distant, has been observed at high resolu- 
tion (A/6A =* 600, 000) by Welty et al. (1994). They observed, 
among 8 total interstellar components, three with velocities be- 
tween -11 and -13.7 km s _1 , with a column density weighted 
average velocity of -12.5 km s _I . The total column density for 
these three components is logA^(Nai) = 12.20 cirT 2 . Because 
7rAqr is more distant than BD +17 4708, several more com- 
ponents are observed, and a higher column density is to be 
expected, but the agreement between the two lines of sight, 
further confirms the interstellar origin of the absorption lines 
observed in the BD +17 4708 spectrum. 

Although we are likely observing several blended ISM 
components toward BD +17 4708, because the absorption is 
optically thin, a single component fit to the entire absorp- 
tion feature should provide a precise measurement of the to- 
tal Nai column density. It has been shown that, even though 



most of the ISM Na is ionized, the total Nai column den- 
sity correlates very well with the total hydrogen column den- 
sity, N(Ui + H 2 ) (Ferlet et al. 1985). Using the relation be- 
tween N(Na i) and N(U I + H 2 ) provided by Ferlet et al. (1985), 
which holds for Nai column densities in the range 10.0 < 
logJV(Nai) cirT 2 < 13.0, we derive a total hydrogen column 
density of logJV(Hi + H 2 ) = 19.78 cirT 2 . Bohlin et al. (1978) 
provide a calibration to transform a total hydrogen column den- 
sity into an E(B - V) value, which for the highest extinctions 
has been confirmed by Rachford et al. (2002). Using this rela- 
tion we obtain a reddening value of E(B - V) — 0.010. 

No formal error bars are provided for the transformations 
from N(Na i) to N(H I + H 2 ) or N(U I + H 2 ) to E(B - V). 
However, we expect the estimate of total hydrogen column 
density from A^(Nai) to be very good because our observed 
Na i column density is right in the middle of the distribution of 
points Ferlet et al. (1985) used to calibrate this relation. The 
transformation from N{H I + H 2 ) to E(B - V) is more difficult 
because we are at the lower end of the distribution of measure- 
ments used by Bohlin et al. (1978). At these low reddenings, 
there is some dispersion in the relation due to the low number 
of discrete absorbers along short sightlines, whereas the rela- 
tion is significantly tighter for more distant lines of sight, over 
which a much larger number of ISM environments are aver- 
aged. The lowest hydrogen column densities used by Bohlin 
et al. (1978) are comparable to the column density we ob- 
serve toward BD +17 4708. Those targets with similar hydro- 
gen column densities have measured reddenings in the range of 
E(B - V) = 0.01 to 0.02. 

Other methods to estimate E(B - V) were also used. The 
E(B - V) value from the Schuster & Nissen (1989) calibration, 
which is based on Stromgren photometry, is essentially zero. 
Interstellar extinction surveys by Fitzgerald (1968) and Arenou 
et al. (1992) suggest E(B - V) = 0.000 and 0.017, respectively, 
when used in conjunction with the Hakkila et al. (1997) code, 
which takes into account the distance to the star. The empirical 
laws by Bond (1980) and Chen et al. (1998) suggest E(B - 
V) = 0.024 while the integrated extinction maps by Burstein 
& Heiles (1978) and Schlegel et al. (1998) set upper limits of 
E(B -V) - 0.043 and 0.035, respectively. Note, however, that 
all the map estimates have large error bars. A simple mean of 
the E(B - V) from the maps results in 0.014 + 0.010. 

Another way to estimate E(B - V) is by means of the 
use of several homogeneously calibrated unreddened color- 
temperature relations. In principle, the standard deviation from 
the mean of several color temperatures minimizes when the 
appropriate E(B - V) value is used. Thus, we used 14 of the 
color-temperature relations by Ramirez & Melendez (2005b) 
and found the dispersion from the mean T e g to be a minimum 
with E(B - V) = 0.008 ±0.001 (see FiggJ. The real error bar in 
this E(B - V) value is very likely to be larger due to photomet- 
ric uncertainties and systematic errors in the r e ff-color calibra- 
tions. Also, the presence of a cool companion (see ^2.51 may be 
affecting this estimate by increasing the red and infrared fluxes 
compared to the case of a single star. In fact, using only blue- 
visible colors we obtain E(B -V) = 0.009 + 0.002 (Fig.gJ. 
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Fig. 3. The observed spectral region around the Na D lines for BD +17 4708, uncorrected for the radial velocity of the star, is 
shown with the histogram (these observational data are described in 33. 1> . The dotted line is the model for telluric water vapor, 
open circles correspond to the model fit to the stellar lines (see 33.51 . and the thick solid line is the model fit to the interstellar 
NaD lines (including blends with telluric lines). 
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Fig. 4. Top panel: mean T e ff obtained from the r e jf-color cal- 
ibrations by Ramirez & Melendez (2005b) as a function of 
E(B - V) using all colors (14, open circles), and blue-visible 
colors only (7, filled circles). Bottom panel: as in the top panel 
for the standard error cr^ = crT e n/ yN, where <x is the standard 
deviation and N the number of colors used. 



In summary, BD +17 4708 is slightly affected by interstel- 
lar reddening but the exact E(B - V) is uncertain. Interstellar 
extinction maps suggest E(B-V) = 0.014±0.010.The£(fi-V) 
value from the r e ff-color relations has a more reasonable er- 
ror of about 0.003, including systematic errors. If we assume 
a slightly smaller error to the E(B - V) value from the fit 
to the ISM lines, e.g., 0.002, we obtain a weighted mean of 
E(B — V) = 0.010 + 0.003. The independent error bars given 
here are somewhat arbitrary but appropriate for the estimate of 
the weighted mean. 

2.2. Effective temperature 

The BL04b data provides a very reliable way to determine the 
effective temperature of BD +17 4708 with the help of theoret- 
ical flux distributions, provided logg, [Fe/H], and E(B - V) 
are known with sufficient accuracy. In particular, the high 



sensitivity of the UV continuum flux to can be used to 
constrain T e $ to a level of about 100 K. Even though the 
completeness of the UV continuum and line opacities in the 
Kurucz models is controversial (Bell et al. 1994, Balachandran 
& Bell 1998), the observed UV fluxes have been shown to 
be reasonably well reproduced by Kurucz models in the Sun 
(Allende Prieto et al. 2003b), Vega (Bohlin & Gilliland 2004a, 
Garcfa-Gil et al. 2005), and late-type stars of different metal- 
licities (Allende Prieto & Lambert 2000). We did not include 
the far UV (A < 2900 A) to quantify the quality of the fits 
in Fig. |2 Not only does this avoid possible errors in the UV 
fluxes of Kurucz models but it also reduces the impact of er- 
rors in the observed flux distribution, which is less accurate 
at short wavelengths. Nonetheless, including the far UV does 
not change our conclusions significantly. We also excluded the 
strong hydrogen lines because they are affected by non-LTE. 

As shown in Fig. |2 the UV and blue-visible spectral re- 
gions are not well fitted with the T e g = 5950 K model. Any 
attempt to reconcile the model and observed UV continuum 
fluxes with an increase in the scaling factor ruins the agree- 
ment in the infrared. The best fit to the data is obtained with 
T eB = 6141 K for E(B - V) = 0.01 (see Fig.|5j». Given the er- 
ror in the extinction value derived in 32.11 the temperature of 
BD +17 4708 is well constrained, from the flux fit, at the 50 K 
level, i.e., T e g - 6141 + 50 K. This 50 K error includes the un- 
certainties in the other atmospheric parameters (logg, [Fe/H]) 
but it is still dominated by the error in the E(B - V) value. A 
systematic error due to the choice of models is certainly present 
but not included in the 50 K. 

Most previous spectroscopic studies of BD +17 4708 have 
found a low r e fr of about 5950 K (see Tableland Fig. ^. 4 
In the Cayrel de Strobel et al. (2001) catalog, for example, 
the 8 entries found for this star have between 5790 K 
and 6100 K, with the mean being 5960 K. Application of 
the InfraRed Flux Method (IRFM) for BD +17 4708 resulted 
in 5955 K according to Ramirez & Melendez (2005a), who 
adopted E(B -V) = 0.000. Note, however, that if E{B -V) = 
0.010 is adopted, the IRFM color-temperature calibrations by 
Ramirez & Melendez (2005b) suggest T e ff 6100 K and if we 



4 By 'spectroscopic' here we refer to T c k obtained from the excita- 
tion and/or ionization balance of iron lines. 
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Fig. 5. Upper panel: mean value of the ratio of observed (/ bs) 
to scaled theoretical (/k m ) fluxes (see Fig. |5J as a function of 
T e ff. Bottom panel: the quality of the fit, as given by x' - 
^X 2 l(n - 1) where y 1 = £?=i[(/obs - /kur)?/of], as a function 
of r e ff (<x is the error in / b s only). 

use blue-visible colors only we obtain T e ff 6150 K (Fig.0J. 
The low temperature obtained from the IRFM and red/infrared 
photometry is probably due to the presence of a cool com- 
panion (see 32. 5> . The three highest temperatures in Table [2 
are those by Axer et al. (1994), Melendez & Ramirez (2004), 
and Asplund et al. (2006). Axer et al. and Asplund et al. de- 
rived their T e g from fitting the wings of the Balmer lines, while 
Melendez & Ramirez used several IRFM temperature-color re- 
lations with a relatively high E(B - V) ^ 0.02 value. 

2.3. Bolometric flux 

The observed absolute flux curve by BL04b covers spectral 
regions that are difficult to model, namely the UV and blue- 
visible. Beyond 1.0 fim, the flux distribution is very well be- 
haved and accurately reproduced by the stellar atmosphere 
models. Furthermore, the infrared portion of the spectrum is 
almost insensitive to the choice of effective temperature when 
the theoretical fluxes are normalized at the Rayleigh-Jeans tail. 

The bolometric flux of BD +17 4708 was obtained by in- 
tegrating the observed flux distribution up to 1.0 fim and the 
predictions from the models for longer wavelengths. Note that 
for the model fits, as in Fig. |2j the theoretical spectra are the 
ones that have been reddened. For the bolometric flux calcula- 
tion, on the other hand, we unreddened the observed flux dis- 
tribution. Thus, we use the term bolometric flux in the intrinsic 
sense, i.e., we use it to refer to the flux that would be measured 
at the top of the Earth's atmosphere in the absence of interstel- 
lar absorption. 

The mean error of each point on the observed flux dis- 
tribution is about 2%. Assuming the error in the models is 
negligible, our best estimate for the bolometric flux is /b i = 
4.89 + 0.10 x 10~ 9 erg crcT 2 s . This value is in good agree- 
ment with that obtained from the Alonso et al. (1995) pho- 
tometric calibrations for the (K, V - K) pair, which result in 
4.80 x 10~ 9 erg cirT 2 s -1 . 



2.4. Angular diameter 

The angular diameter of BD +17 4708 can be calculated with 
the /boi and derived above. This value corresponds to the 
limb-darkened angular diameter. Propagating the 2% error in 
/boi and 1 % error in T e $ we determine our best solution for the 
angular diameter as: = 0.1016 + 0.0023 mas. 

2.5. Binarity 

Our model fits (and hence our results for /boh T e ff, etc.) con- 
sider BD +17 4708 as a single star. However, it is known 
that BD +17 4708 shows periodic radial velocity variations 
with an amplitude of about 4.7 km s _1 over a period of about 
220 days (Latham et al. 1988). Model fits to the radial veloc- 
ity curve by Latham et al. suggest a mass function f(M) - 
(m 2 sin 3 /(wi + m 2 ) 2 = 0.0019 + 0.0004M Q , where m,\ and m 2 
are the masses of the stars and i the orbital inclination. 

A rough estimate assuming a mass of about 0.9M© for the 
primary (see 33.5. Il l, results in a companion mass of O.15M 
(using < sin 3 i >- 3/5), which corresponds to a late M type 
star with a T e s ~ 3000 K. Using Kurucz model fluxes we find 
that the contribution of the secondary to the UV and blue fluxes, 
which are much more sensitive to T e ff than the IR, is negligible 
(less than 1/50). Although the secondary contributes a signif- 
icant flux in the IR (about 1/5), the shape of the spectral en- 
ergy distribution is nearly unchanged. Given that it is this shape 
along with a scaling factor what determine the best r e ff solu- 
tion for the primary, including the companion flux in the fits 
will not affect the r e ff result significantly. If included, the scal- 
ing factor would need to be reduced and the r e ff of the primary 
increased to match the observed UV and blue fluxes. However, 
given that the orbital inclination is unknown, it is safer to use 
a single model flux to fit the observed energy distribution, but 
note that the companion flux may have an important effect in 
the observed colors, making them redder. This is probably the 
reason why direct application of the IRFM suggest a lower T e ff 
for the primary compared to the r e ff obtained from the flux fit 
(£3- 

The IRFM temperature of this star is 5950 K according to 
Ramirez & Melendez (2005a), who used E(B - V) = 0. Using 
E(B - V) = 0.01 the IRFM temperature increases to 6025 K, 
about 120 K lower than the T e ff obtained from the flux fit. It is 
unlikely that such large difference is due to errors in the abso- 
lute infrared flux calibration and/or the zero point determina- 
tion of the IRFM scale (see Ramirez & Melendez 2005a 
for details). The most likely reason for this discrepancy is the 
flux contributed by the companion, which is more important in 
the infrared. If the flux at a given wavelength in the infrared 
(/ir) is larger then the ratio R = /boi//iR is smaller compared 
to that for a single star. This /^-factor is the T e s indicator in 
the IRFM, roughly proportional to (Ramirez & Melendez 
2005a). Thus, the IRFM temperature obtained for a star with an 
ignored cool companion is underestimated. In order to account 
for the 120 K difference (about 2% error in T e ff), an error of 
about 6% in the 7?-factor is required. In the previous paragraph 
we estimated a 20% extra infrared flux due to the companion. 
Using the same models, the bolometric flux increases by about 
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10% if the cool companion is included. If this is the case, the 
/^-factor has been underestimated by about 8%. However, note 
again that these flux estimates are not accurate due to the large 
uncertainty in the mass and temperature of the companion. 

3. Spectral line analysis 

3.1. Observations 

BD +17 4708 was observed from McDonald Observatory on 
October 30, 2004 UT using the 2dcoude spectrograph (Tull 
et al. 1995) and the 2.7m Harlan J. Smith telescope. Four in- 
dividual exposures of 20 min each were obtained at the fo- 
cal station F3 using grating E2 -a 53.67 gr mrrr 1 R2 echelle 
from Milton Roy Co.-, a 1.2 arcsec slit, and a 2048 x 2048 
Tektronix CCD. The spectra have a FWHM resolving power 
of A/ 6 A ^ 60, 000 with full spectral coverage from 3600 A to 
5300 A, and substantial but incomplete coverage from 5300 A 
to 10000 A. The spectra were reduced using the echelle pack- 
age in IRAF. 5 The bias level in the overscan area was modeled 
with a polynomial and subtracted. An ultra-high signal-to-noise 
flatfield was used to correct pixel-to-pixel sensitivity variations, 
and the scattered light was modeled with smooth functions and 
removed. 

The spectra were optimally extracted after cosmic-ray 
cleanup, and calibrated in wavelength with a Th-Ar hollow 
cathode lamp (Allende Prieto 2001). By cross-correlating the 
four individual spectra, we concluded that shifts among them 
were smaller than 0.2 km s _1 , and we simply coadded their 
signal obtaining a single spectrum with a signal-to-noise ratio 
per pixel in excess of 300 between 5000 A and 8000 A, and in 
excess of 100 between 4000 A and 10000 A. The individual or- 
ders were continuum normalized, combining the signal for the 
wavelength intervals registered in multiple orders. Similarly to 
the procedure described by Barklem et al. (2002), we took ad- 
vantage of the slow variation of the blaze function between or- 
ders in the normalization process, in order to derive reliable 
line shapes for the strongest lines. 

3.2. Atomic data 
3.2.1. Iron 

All the gf values for the iron lines used in this work have been 
measured in the laboratory. No attempts to reduce the line-to- 
line scatter in the abundances using differential analysis or as- 
trophysical gf values have been made. Thus, our derived iron 
abundances are strictly given on an absolute scale. 

The original sources for the transition probabilities of the 
Fei lines are listed by Lambert et al. (1996), who extensively 
compared them and concluded that they were all essentially on 
the same scale, although minor corrections are needed in a few 
cases. The gf values for the Fe n lines have been adopted from 

5 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation - http://iraf.noao.edu 



Melendez et al. (2006), who use gf values from theoretical cal- 
culations put onto the laboratory scale by means of laboratory 
lifetimes and branching ratios. Note that these ^/-values are 
very similar to those in Lambert et al. (1996). The mean dif- 
ference in log gf, in the sense Melendez et al. 2006 - Lambert 
et al. 1996, for 4 Fe n lines available in both studies and given in 
Table|2] is only 0.03 dex. A similar comparison with the com- 
pilation by Allende Prieto et al. (2002) reveals that, on average, 
their gf values are on the same scale as those by Melendez 
et al. (2006). However, the line-to-line scatter reduces when 
adopting the latter set of gf values. 

Regarding van der Waals pressure broadening, almost all 
the damping constants adopted in this work are from Barklem 
et al. (2000) and Barklem & Aspelund-Johansson (2005). For 
a few Fe i lines not included in the Barklem et al. tables, the 
classical Unsold approximation, enhanced by a factor of 2, was 
adopted. Standard radiative (e.g., Gray 1992) and Stark (e.g., 
Cowley 1971) broadening approximations, as coded in the lat- 
est version of MOOG (Sneden 1973), 6 were used. 

The iron line data and equivalent widths (E.W.s) measured 
in the spectrum of BD +17 4708 are given in Tabled Gaussian 
profile fitting was used to measure the line E.W.s. 

3.2.2. Other elements 

For the strong 5180 A Mgib lines, as well as for the 7774 A 
O i triplet, transition probabilities were obtained from the NIST 
Atomic Spectra Database. 7 

It was difficult to find weak Mg lines with reliable transition 
probabilities in our spectrum. In fact, we found only the un- 
blended weak 4571 A Mgi line, with a reliable gf value from 
the NIST database. Another weak Mg i line is that at 571 1 A, 
for which we used the solar gf derived by Fuhrmann et al. 
(1995). 

Data for four weak Si i lines were taken from the compi- 
lation by Allende Prieto et al. (2004a), who concentrated on 
lines with transition probabilities measured in the laboratory or 
obtained from accurate theoretical calculations. 

For the Ca abundance determination, we used the line list 
by Bensby et al. (2003) but adopting the gf values from the 
NIST database instead of using their solar g/'s. The Bensby 
et al. gf values are systematically lower by about 0.2 dex com- 
pared to those obtained from the NIST database. However, the 
line-to-line scatter in our derived mean Ca abundance is similar 
for the two sets of gf values. 

Radiative, Stark, and van der Waals broadening was com- 
puted in the same way as for the iron lines. Note that in this case 
all lines are weak (with the exception of the Mg i b triplet and 
the 6439 A Ca i line) so the use of the modified Unsold approx- 
imation, when necessary, to obtain the van der Walls damping 
constants instead of using those from the theory of Barklem 
et al. (2000) has no noticeable effect on the abundances. 



6 http ://verdi .as.utexas. edu/moog . html 

7 http://physics.nist.gov/PhysRefData/ASD/lines_form.html 
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Wavelength Species E.R logg/ logT (1— or)/2 E.W. 

A eV md cm 3 s -1 mA 



Table 2. Iron line data. T and a are the Van der Waals FWHM 
per perturber at 10,000 K and velocity parameter, respectively. 
The damping constants for the lines marked with a t corre- 
spond to the modified Unsold approximation. The last col- 
umn gives the equivalent widths measured in the spectrum of 
BD +17 4708. 



3.3. Modeling 

Spectrum synthesis was performed using MOOG (Sneden 
1973) and the non-LTE codes TLUSTY and SYNSPEC 
(Hubeny 1988, Hubeny & Lanz 1995). 8 For practical reasons, 

8 http://nova.astro.umd.edu 



Wavelength Species E.R logg/ logT (l-cc)/2 E.W. 

A eV rad cm 3 s~' mA 



7771.944 


Oi 


9.146 


0.37 


-7.469 


0.383 


38.2 


7774.166 


Oi 


9.146 


0.22 


-7.469 


0.383 


31.7 


7775.388 


Oi 


9.146 


0.00 


-7.469 


0.383 


22.9 


4571.096 


Mgi 


0.000 


-5.39 


-7.645 


0.377 


19.4 


5167.321 


Mgi 


2.709 


-0.86 


-7.267 


0.381 





5172.684 


Mgt 


2.712 


-0.38 


-7.267 


0.381 





5183.604 


Mgi 


2.717 


-0.16 


-7.267 


0.381 





5711.100 


Mgi 


4.346 


-1.67 


-7.218 


0.300 t 


17.6 


5708.397 


Sii 


4.954 


-1.37 


-7.183 


0.300 t 


10.0 


5948.540 


Si i 


5.082 


-1.13 


-7.169 


0.300 t 


12.0 


7918.382 


Sii 


5.954 


-0.51 


-7.010 


0.300 t 


10.3 


7932.348 


Si i 


5.964 


-0.37 


-7.006 


0.300 t 


10.5 


4526.928 


Cat 


2.709 


-0.42 


-7.021 


0.300 t 


15.8 


4578.551 


Cat 


2.521 


-0.56 


-7.125 


0.300 t 


18.3 


5512.980 


Cat 


2.933 


-0.30 


-7.269 


0.300 t 


16.7 


6166.439 


Cai 


2.521 


-0.90 


-7.146 


0.372 


1 1.3 


6169.042 


Cat 


2.523 


-0.54 


-7.146 


0.372 


21.0 


6169.563 


Cat 


2.526 


-0.27 


-7.145 


0.372 


30.6 


6439.075 


Cat 


2.526 


0.47 


-7.569 


0.379 


69.4 


6471.662 


Cat 


2.526 


-0.59 


-7.570 


0.380 


24.4 


6493.781 


Cat 


2.521 


0.14 


-7.571 


0.381 


54.8 


6499.650 


Cai 


2.523 


-0.59 


-7.571 


0.381 


18.1 



Table 3. As in Table for the O, Mg, Si, and Ca fines. 
Equivalent widths are given for all but the strong Mg ib lines. 



MOOG was preferred for matching the equivalent widths of the 
iron lines, while SYNSPEC was used to fit the profiles of strong 
lines. The same scaled solar abundances (those by Grevesse & 
Sauval 1998) were used in the two codes and thus, only very 
small differences, mainly due to the continuum opacity calcu- 
lations, may be present when comparing the results from the 
two codes. All line formation calculations were done assuming 
LTE, with the only exception of the 7774 A triplet (see |4j. 

For the line-profile fitting, all the synthetic profiles have 
been broadened by convolving the theoretical spectra with 
Gaussian profiles of FWHM = 0.18 A in the red and FWHM = 
0.21 A in the near infrared. These FWHM values are empir- 
ically determined global broadening parameters that fit very 
well weak lines with reliable atomic data. In fact, given the 
resolution R ~ 60,000 (6 A ~ 0.09 A at 5500 A, 6 A. ~ 0.13 A 
at 7780 A), a solar-like macroturbulent velocity of 1.5 km s _1 
(5A ~ 0.03 A at 5500 A, 5A ~ 0.04 A at 7780 A), and a low 
projected rotational velocity of v sin i ~ 3 km s , our estimates 
for the FWHM values are well justified. 

We used the most recent Kurucz no-overshoot model at- 
mospheres with of-element enhancement (e.g., Kurucz 1970, 
1979). The use of models with the convective overshooting op- 
tion switched on produces an almost constant shift of less than 
0. 1 dex in the abundance scale but preserves abundance ratios 
as well as the difference in the mean Fe abundances from Fe i 
and Fe n lines reported in ^3.61 Use of models with solar scaled 
abundances (i.e., without c-element enhancement) produced 
essentially the same abundances. For each set of atmospheric 
parameters adopted, a microturbulent velocity v t was derived 
by making the abundances from the Fe i lines independent of 
their reduced equivalent widths (E.W. I A). 

3.4. T e ff from the Balmer lines 

Balmer line-profiles were synthesized using the prescription by 
Barklem et al. (2002) but adopting Kurucz model atmospheres. 
In short, Stark broadening was computed according to Stehle & 



4630.120 


Fei 


2.279 


-2.52 


-7.518 


0.373 


7.0 


4745.800 


Fei 


3.654 


-1.27 


-7.356 


0.300 t 


8.6 


4939.686 


Fei 


0.859 


-3.34 


-7.748 


0.377 


23.9 


4994.129 


Fei 


0.915 


-3.07 


-7.744 


0.377 


28.0 


5012.068 


Fei 


0.859 


-2.64 


-7.751 


0.377 


53.0 


5044.211 


Fei 


2.851 


-2.03 


-7.280 


0.381 


8.3 


5051.634 


Fei 


0.915 


-2.79 


-7.746 


0.377 


42.3 


5079.740 


Fei 


0.990 


-3.22 


-7.739 


0.378 


20.1 


5110.413 


Fei 


0.000 


-3.76 


-7.826 


0.373 


42.4 


5123.720 


Fei 


1.011 


-3.07 


-7.739 


0.378 


26.0 


5127.359 


Fei 


0.915 


-3.31 


-7.749 


0.377 


19.9 


5150.839 


Fei 


0.990 


-3.00 


-7.742 


0.377 


23.0 


5151.911 


Fei 


1.011 


-3.32 


-7.740 


0.377 


17.3 


5166.282 


Fei 


0.000 


-4.20 


-7.827 


0.373 


21.2 


5171.596 


Fei 


1.485 


-1.78 


-7.688 


0.373 


63.4 


5194.941 


Fei 


1.557 


-2.08 


-7.680 


0.373 


45.1 


5216.273 


Fei 


1.608 


-2.14 


-7.674 


0.372 


39.4 


5227.189 


Fei 


1.557 


-1.23 


-7.681 


0.373 


90. 1 


5228.376 


Fei 


4.220 


-1.19 


-7.233 


0.361 


4.4 


5253.461 


Fei 


3.283 


-1.57 


-7.203 


0.386 


6.2 


5307.360 


Fei 


1.608 


-2.98 


-7.678 


0.373 


1 1.8 


5322.041 


Fei 


2.279 


-2.89 


-7.600 


0.382 


4.2 


5328.531 


Fei 


1.557 


-1.85 


-7.685 


0.374 


62.0 


5332.899 


Fei 


1.557 


-2.78 


-7.685 


0.374 


16.0 


5341.023 


Fei 


1.608 


-1.95 


-7.679 


0.373 


52.5 


5371.489 


Fei 


0.958 


-1.65 


-7.753 


0.376 


95.9 


5373.708 


Fei 


4.473 


-0.74 


-7.123 


0.359 


6.5 


5397.127 


Fei 


0.915 


-1.99 


-7.759 


0.375 


81.7 


5429.696 


Fei 


0.958 


-1.88 


-7.755 


0.376 


87.8 


5432.947 


Fei 


4.445 


-0.94 


-7.153 


0.360 


7.2 


5434.523 


Fei 


1.011 


-2.12 


-7.750 


0.377 


72.5 


5473.900 


Fei 


4.154 


-0.72 


-7.266 


0.380 


8.4 


5497.516 


Fei 


1.011 


-2.85 


-7.752 


0.376 


41.0 


5501.465 


Fei 


0.958 


-3.04 


-7.757 


0.375 


32.4 


5506.779 


Fei 


0.990 


-2.80 


-7.754 


0.376 


42.1 


5543.935 


Fei 


4.217 


-1.04 


-7.263 


0.381 


5.1 


5638.262 


Fei 


4.220 


-0.77 


-7.270 


0.382 


6.8 


5701.544 


Fei 


2.559 


-2.22 


-7.576 


0.382 


8.5 


5905.671 


Fei 


4.652 


-0.69 


-7.144 


0.359 


6.0 


5930.179 


Fei 


4.652 


-0.17 


-7.149 


0.359 


24.8 


5934.654 


Fei 


3.928 


-1.07 


-7.153 


0.377 


7.6 


6003.012 


Fei 


3.881 


-1.06 


-7.181 


0.380 


7.8 


6027.050 


Fei 


4.076 


-1.09 


-7.397 


0.300 t 


6.3 


6056.004 


Fei 


4.733 


-0.40 


-7.130 


0.357 


4.7 


6170.507 


Fei 


4.795 


-0.38 


-7.119 


0.355 


6.0 


6200.313 


Fei 


2.608 


-2.44 


-7.588 


0.382 


5.0 


6213.430 


Fei 


2.223 


-2.48 


-7.691 


0.368 


8.9 


6232.641 


Fei 


3.654 


-1.22 


-7.498 


0.300 f 


7.9 


6265.133 


Fei 


2.176 


-2.55 


-7.699 


0.369 


1 1.2 


6344.148 


Fei 


2.433 


-2.92 


-7.620 


0.377 


3.9 


6419.949 


Fei 


4.733 


-0.24 


-7.193 


0.363 


12.0 


6609.110 


Fei 


2.559 


-2.69 


-7.610 


0.377 


4.2 


6750.152 


Fei 


2.424 


-2.62 


-7.608 


0.380 


6.9 


6841.338 


Fei 


4.607 


-0.71 


-7.258 


0.367 


5.7 


6855.162 


Fei 


4.558 


-0.74 


-7.347 


0.300 t 


5.4 


7090.383 


Fei 


4.230 


-1.11 


-7.165 


0.376 


5.0 


4620.521 


Fe n 


2.828 


-3.21 


-7.878 


0.347 


7.4 


4629.339 


Fe ii 


2.807 


-2.28 


-7.886 


0.372 


42.0 


5197.577 


Fe ii 


3.230 


-2.22 


-7.881 


0.377 


32.5 


5234.625 


Fe ii 


3.221 


-2.18 


-7.881 


0.376 


38.0 


5264.812 


Fe ii 


3.230 


-3.13 


-7.875 


0.350 


7.1 


6432.680 


Fe ii 


2.891 


-3.57 


-7.899 


0.398 


5.6 


6516.081 


Fe ii 


2.891 


-3.31 


-7.899 


0.399 


8.8 
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Fig. 6. Top panel: the observed Ha profile (points) 
is shown along with model line profiles of r e ff = 
5950,6050,6150,6250,6350,6450 K (solid lines). The 
bars on top show the regions excluded from the x 2 calculation. 
Bottom panel: quality of the fits in the top panel, as measured 

by*' = y/x 2 /(n - 1) where x 2 = [(/<*» - f^) 2 /ojl (a is 
the error in / Q b s only), as a function of T e g (filled circles). The 
solid line is a cubic fit to the filled circles. 

Hutcheon (1999) while self-broadening is from Barklem et al. 
(2000). The shapes of strong lines like Ha are very well deter- 
mined in our spectrum by fitting the blaze shape of each clean 
order (those free from very strong line absorption), and mod- 
eling the smooth variation of the shape of the blaze with order 
number to set the continuum (see §2 in Barklem et al. 2002 for 
details). 

The results for Ha are shown in Fig. |6] along with the 
observed profile in the spectrum of BD +17 4708. As shown 
in Fig. [6] the wings of the Ha line are very sensitive to r e ff. 
Unfortunately, due to contamination by metallic lines, the H/3 
profile is difficult to use as a temperature indicator. 

Excluding the most prominent metallic features, as shown 
in Fig. U test favors a r e ff = 6165 K from Ha. For the 
X 1 calculations, the Ha line was cut at 6550 A since shorter 
wavelengths fall outside the CCD. In fact, our spectrum in this 
order goes down to about 6535 A, but the data points between 
6535 A and 6550 A are, observationally, somewhat more un- 
certain. If the blue wing of the Ha line is ignored altogether (to 
minimize further the observational errors), the temperature in- 
creases only by 20 K. On the other hand, the sensitivity to T e g 
decreases for A > 6590 A and, given the S/N, introducing these 
longer wavelengths to assess the quality of the fits increases 
only the absolute x 2 values without changing significantly the 
inferred T e g. 

Our Ha temperature is in excellent agreement with that 
given by Asplund et al. (2006), who obtained T e g = 6183 K 



from fits to their Ha profile using MARCS models but the same 
treatment for the line broadening. Both, ours and Asplund et al. 
Ha temperatures are in good agreement with the T e g derived 
from the fitting of the spectral energy distribution. 

3.5. Surface gravity from strong lines 

The stellar surface gravity of a star can, in principle, be 
obtained from an estimate of its mass and its measured 
trigonometric parallax, besides reasonable estimates of T e s and 
[Fe/H]. The mass of a nearby star can be reasonably estimated 
from its position on a color-magnitude diagram using theo- 
retical isochrones but the Hipparcos parallax of stars farther 
than 100 pc, as is the case for BD +17 4708, is quite uncertain 
and therefore their trigonometric logg values are not reliable. 
In fact, using this method we only obtain a weak constrain: 
3.8 < logg < 4.6. 

Fortunately, the wings of some strong lines are sensitive to 
the log g value and are less affected by T e g and can thus be used 
to constrain the surface gravity. In the BD +17 4708 spectrum, 
only the two strongest lines of the Mg i b triplet at 5 172.7 A and 
5183.6 A seem suitable for this kind of analysis (Fig.0. Note 
that the cores of the strong lines are strongly affected by non- 
LTE and we do not expect good fits in the line centers, only the 
wings should be used to assess the quality of the fits. 

Fig.0shows that an excellent fit to the wings of the Mg ib 
triplet is obtained with logg = 3.87 when the Mg abundance 
is set to AMg = 6.19, as derived in ^3.71 Our preferred solution 
for the surface gravity is thus logg = 3.87 ± 0.08. The error 
bar was estimated by propagating the error of 0.06 dex in A Mg , 
which includes the 50 K error in T e g . 

The Mg i b triplet is very strong and contaminated by metal- 
lic lines in the solar spectrum. Nevertheless, adopting the same 
procedure we used to obtain the \ogg value of BD +17 4708, 
we were able to satisfactorily reproduce the wings of these lines 
in the solar spectrum of Kurucz et al. (1984) with the standard 
solar logg = 4.44 and Mg abundance of Af^g = 7.53 (Asplund 
et al. 2005). A visual inspection showed that the accuracy of 
our method of log g determination in the Sun is about 0. 1 dex. 
Therefore, the \ogg value we derive for BD +17 4708 is accu- 
rate at the 0. 1 dex level in the absolute scale. 

3.5.1 . Mass, age, and radius from theoretical 
isochrones 

Isochrones in the theoretical HR diagram (T e g vs. log g) instead 
of the observational HR diagram (absolute magnitude vs. color) 
can be used to estimate the mass and age of a star if its parame- 
ters, but not necessarily its distance, are known with accuracy. 
This is the case of BD +17 4708. 

Although the mass estimates are normally accurate using 
this approach, the age determinations may be subject to se- 
vere systematic errors and statistical biases (see, e.g., Pont & 
Eyer 2004) so they must not be considered accurate even if 
the stellar parameters are. We used the Bertelli et al. (2004) 
isochrones, as in Allende Prieto et al. (2004a), to estimate the 
mass (M) and age (f) of BD +17 4708. The Bertelli et al. 
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5183.0 5183.5 5184.0 5184.5 
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Fig. 7. Top panel: the region around the strong Mglb triplet as given by the observations (open circles), and as predicted by 
models of T e g = 6141 K, [Fe/H] = -1.74, Am s = 6.19, and three different values of the surface gravity (logg = 3.5, 3.87,4.5). 
Residuals of the fit are also shown. Bottom panel: zoom into the 5183.6 A line. 



isochrones were computed using solar-scaled chemical compo- 
sitions. In metal-poor stars, however, the or-element enhance- 
ment ([of/Fe] ^ +0.4 in our case, 33. 7> has an important ef- 
fect on these calculations (e.g., VandenBerg et al. 2000, Kim 
et al. 2002) although its effect on the mass and age derived 
from the isochrones is relatively small (about +0.1 M and 
-0.4 Gyr in our case). Therefore, we increased the [Fe/H] 
value of BD +17 4708 by about 0.2 dex to mimic the a-element 
enhancement, as suggested by Salaris et al. (1993) and obtained 
M = 0.91+°;i'M Q and t = 8.8;^ Gyr (2cr errors). 

The referee noted that using our derived parameters the 
Vandenberg et al. (2000) isochrones suggest an age close to 10 
Gyr but if the log g value is increased to 4.05 then the age would 
increase to about 13.5 Gyr. Using the Bertelli et al. isochrones 
and logg = 4.05, we obtained t = 11.4^2 Gyr, i.e., an in- 
crease of 2.6 Gyr in the mean age. The halo is believed to have 
an age of about 13 Gyr (e.g., Schuster et al. 2006). However, 
given that this mean age is calculated using large samples of 
halo stars and in some cases sophisticated statistics, this should 
not be used to discard or confirm ages of individual stars. Note 
that, for example, our derived age is in agreement with that 
given by Nordstrom et al. (2004), who took into account the 
statistical biases in isochrone age determinations described in 
Pont & Eyer (2004). Also, some halo stars, even more metal- 
poor than BD +17 4708, seem to be younger than the mean age 
of the halo (see, e.g., Table 2 in Li & Zhao 2004, who give a 
compilation of radioactive ages, including theirs). 

The radius that we obtain using isochrones is about 1 .8/? Q , 
with a 2-cr range that goes from 1 .5 to 23R if we include sys- 
tematic errors in our T e s and logg estimates. Although inaccu- 
rate, the Hipparcos parallax constrains the radius to a 1-cr range 
from 1.1 to 1.6 R G , if we adopt our angular diameter ( 32.4> . A 



slightly higher log g value, for example log g — 4.0, would re- 
sult in R ^ 1 AR Q . Note that this would still be in good, albeit 
marginal, agreement with our result for log g considering the 
random error bar (0.08 dex) and a possible systematic error in 
the absolute scale (about 0.10 dex). 

3.6. The iron abundance 

The customary approach to determination of the Fe abundance 
invokes LTE for the excitation and ionization of iron neutral 
atoms and singly-charged ions. An estimate of the effective 
temperature is obtained by the requirement that the derived Fe 
abundances from the Fe i lines be independent of their exci- 
tation potential. Application of this requirement generally de- 
mands a prior determination of the microturbulence (v,), often 
from the same set of Fe i lines and the condition that the Fe 
abundance be independent of a line's E.W. Then, the imposi- 
tion of ionization equilibrium through the requirement that the 
Fe i and Fe n lines give the same Fe abundance defines a locus 
in the (T e ff, log g) plane which with the T e ff from the Fe i lines 
(or another source) serves to determine the surface gravity. 

We used 56 Fe i lines covering the excitation potential (E.R) 
range from to 5 eV and 7 Fe n lines to derive the iron abun- 
dance for various choices of atmospheric parameters (Figs. [S] 
and|5}. The lines we selected have E.W. between 4 and 100 mA 
to avoid errors due to noise and saturation. 

The Fe i lines with an almost 5 e V range in excitation poten- 
tial demand an effective temperature only slightly larger than 
about 6141 K (Fig. [8}, almost independently of the choice of 
surface gravity, and confirming the temperature provided by 
the flux distribution. A temperature of 5950 K is demonstrably 
too low. 
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Fig. 8. Abundance of iron from Fe 1 (open circles) and Fe n lines (filled circles), as determined using different combinations of 
r e ff, log g, and v f ; as a function of excitation potential and reduced equivalent width. The solid lines are linear fits to the Fe 1 data 
only. 



A correlation between the excitation potential of the lines 
and their reduced equivalent widths may lead to degenerate so- 
lutions for the (T^ , v,) pair. There is no such correlation for our 
Fe 1 lines with E.P. > 2 eV. Our Fe 1 lines with E.P. < 2 eV do 
not show such correlation either but they are all stronger (i.e., 
they all have larger reduced equivalent widths) than those with 
E.P. > 2 eV. The strong, low E.P. lines allow to better deter- 
mine v,. Notice, however, that the A Fe vs. E.P. relations shown 
in Fig.|S]do not change dramatically if the lowest E.P. lines are 
avoided. 

The LTE ionization equilibrium is satisfied at the locus 
shown in Fig.^| At T e s = 6141 K, the locus (see also Fig.|9jl 
corresponds to about \ogg = 4.3, a value higher than that pro- 
vided by the fit to the Mg 1 b lines. Note that this result is incon- 
sistent with that found by Edvardsson (1988), who concludes 
that the strong line gravities are larger than those obtained from 
the ionization balance in a sample of subgiants with metallic- 
ities higher than about [Fe/H] = -0.5. At that lower value of 
logg = 3.87 and T e s = 6141 K, the Fe abundance from Fen 
lines is about 0.15 dex less than that from the Fei lines, as it 
is clearly seen in the middle panels of Fig.[S] Ionization equi- 
librium at logg = 3.87 is achieved if T e g 5900 K, but this 
temperature is judged to be too low. 

Our analysis shows that a consistent analysis of the flux 
distribution and the Mg 1 b, Fe 1, and Fe 11 lines cannot be found 
within the constraints of a classical LTE model atmosphere 
analysis. The inconsistencies, almost certainly, cannot be as- 
cribed to the accumulation of errors in the flux and line data. 
One must suspect a failure of the classical atmosphere and/or 
the breakdown of the LTE assumption. 



Introduction of departures from LTE into the formation of 
iron lines within a classical model atmosphere constructed as- 
suming LTE for all sources of continuous and line opacity calls 
for atomic data on radiative and collisional processes far be- 
yond the restricted need for the corresponding LTE analysis. 
The main non-LTE effect on the Fe lines has been shown to 
be an overionization of neutral Fe atoms resulting from the 
UV flux (e.g., Athay & Lites 1972). This effect is the more 
severe for metal-poor stars owing, principally, to the reduced 
line blocking in the UV. 

Calculations reported for HD 140283, a star more metal- 
poor and cooler than BD +17 4708, show that, using ID model 
atmospheres (see next paragraph), the Fe abundance from the 
Fe 1 lines might be increased by up to about 0.5 dex for a non- 
LTE analysis (e.g., Korn et al. 2003, Shchukina et al. 2005) 
while leaving the Fe abundance from the Fe n lines nearly un- 
changed. If these non-LTE effects are taken into account in our 
case, they would increase further the difference between the 
Fe abundance from the neutral and ionized lines and would re- 
quire even higher surface gravities to achieve ionization equi- 
librium. Note, however, that the role of inelastic collisions with 
neutral hydrogen in the non-LTE calculations needs to be ex- 
plored in more detail given that they may significantly reduce 
the size of the non-LTE corrections to the Fe 1 abundance (Korn 
et al. 2003). In fact, we find that they are very important for the 
oxygen abundance determination from the 7774 A triplet al- 
though the formulation commonly adopted (that by Steenbock 
& Holweger 1984) is questionable and the resulting abun- 
dances uncertain (see 

Classical atmospheres with their assumption of plane paral- 
lel homogeneous layers in hydrostatic equilibrium cannot rep- 
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resent the stellar granulation resulting from convective instabil- 
ities. Modeling of granulation in main-sequence stars including 
metal-poor examples is beginning with attendant analyses of 
LTE and non-LTE line formation (Asplund 2005). The latter 
models are commonly referred to as 3D models with classical 
models as ID models. Calculations for HD 140283 suggest that 
the non-LTE Fe abundance from the Fe n lines is increased by 
about 0.3 dex in going from a ID to a 3D model of the same 
atmospheric parameters but the non-LTE abundance from Fe i 
lines is unchanged (Shchukina et al. 2005). 

Taking the results for HD 140283 at face value, and apply- 
ing them to our case, the switch from the ID LTE analysis to 
the 3D non-LTE one (without H collisions) would lead to an 
increase of about +0.2 dex in the Fe I - Fe II abundance dif- 
ference (the Fei abundance increasing by 0.5 dex due to non- 
LTE in ID and the Fen abundance increasing by 0.3 dex due 
to 3D effects in non-LTE). Thus, if the sense and approximate 
magnitude of these effects applies to BD +17 4708, the dis- 
agreement between the surface gravity derived from the LTE 
analysis of the Mg i b lines and Fe ionization equilibrium is in- 
creased. Excitation by collisions with H atoms could alleviate 
the disagreement discussed somewhat but no reliable theory to 
include them in the non-LTE calculations is available at present 
(see Q. 

There is clearly a need for a fuller exploration of the non- 
LTE effects (in particular H collisions) both in the construction 
of the model atmosphere and in the line formation. Detailed 
testing of 3D models is a necessity with confrontation between 
predictions of the energy distribution, line strengths, wave- 
lengths, and asymmetries. Pending this major challenge, we 
conclude that the best fits to the data so far have been achieved 
with r e ff = 6141 K and logg = 3.87. Given that the Fe n lines 
seem to be less affected by errors in or non-LTE (see, e.g., 
Thevenin & Idiart 1999, but see also Shchukina et al. 2005), we 
disregard the ionization balance condition and adopt the mean 
abundance from the Fen lines only as our [Fe/H] indicator. 
Thus, our preferred solution for the metallicity of BD +17 4708 
is [Fe/H] = -1.74 + 0.09, which, for our inferred solar iron 
abundance, 9 corresponds to Ap e - 5.11 ± 0.09. 10 

3. 7. Mg, Si, and Ca abundances 

Besides their importance in the study of Galactic chemical evo- 
lution, the abundances of these elements are relevant for stel- 
lar structure and evolution calculations (e.g., VandenBerg et al. 
2000). Although they are less abundant than carbon, nitrogen, 
and oxygen, they are important sources of opacity at high tem- 
peratures. 

It is well established that the [a/Fe] ratios, where a is an a 
element (e.g., Mg, Si, Ca, etc.), in the majority of metal-poor 
stars are well above solar. Most authors agree on [a/Fe] ratios 

9 Our line list and adopted atomic data result in Af c ,o = 7.51 ± 0.08, 
with no significant difference between the Fei and Fen abundances, 
although not all the lines could be used in the solar spectrum due to 
saturation. Details will be given in Ramirez et al. (2006). 

10 Unless otherwise noted, all the error bars for the abundances are 
l-cr errors. 
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Fig. 9. Mean abundance of iron from Fe i (open symbols) and 
Fe ii (filled symbols) lines as a function of T e s and log g. Here 
[Fe/H] = Ap e - 7.51, i.e., we used A Fe = 7.51 for the Sun. 
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Fig. 10. Locus of the T e ff and log g values for which ionization 
equilibrium of Fe lines is satisfied. The circle with the error 
bars correspond to our best estimates of r e fr from the flux dis- 
tribution and log g from the fit of the Mg i b lines. 



of about +0.3 + 0.1 dex in the halo (e.g., Carretta et al. 2000, 
Gratton et al. 2000, Idiart & Thevenin 2000, Arnone et al. 2005, 
Barklem et al. 2005). 

The Mg abundance was obtained by averaging the abun- 
dance derived from the 4571 A and 5711 A lines. However, 
due to the more reliable transition probability of the former, we 
gave a double weight to the Mg abundance obtained from the 
4571A line. In this way we find A Mg = 6.19 ± 0.05. The four 
Si lines listed in Table[3]resulted in As; = 6.12 ± 0.06 while the 
Ca lines suggest A Ca = 4.93 ± 0.06. 

Almost all the Mg, Si, and Ca lines used to derive the 
abundances given above are very strong (E.W. >~ 90 mA) 
in the solar spectrum, which makes them unsuitable to derive 
solar abundances due to saturation effects. Instead, we used 
the solar abundances derived by Asplund et al. (2005) to ob- 
tain [Mg/Fe] = 0.40 + 0.10, [Si/Fe] = 0.35 ± 0.11, and 
[Ca/Fe] = 0.36 + 0.11, where the error due to the uncertain- 
ties in [Fe/H] and r e ff have been included. The mean [a/Fe] 
ratio is 0.37 + 0.06, where the error bar here is a standard error. 
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Fig. 11. Observed profile of the Oi triplet in BD +17 4708 
(filled circles). Non-LTE model fits to the data (solid lines) are 
shown for S h = 0,1, and 10. LTE fits are also shown. The 
abundances adopted in each case, given in each panel, have 
been chosen to fit the reddest line of the triplet. Residuals, 
shifted by +0.80, are show at the bottom of each panel. 



4. The Oi triplet: non-LTE effects and [O/Fe] ratio 

The observed oxygen abundances are relevant for studies of 
chemical evolution of the Galaxy and supernovae yields (e.g., 
Wheeler et al. 1989), as well as for the modeling of stellar 
structure and evolution (e.g., VandenBerg & Stetson 1991). 

The IR triplet lines are strong enough as to be detected in 
most metal-poor FGK dwarfs and it has been long known that 
they suffer from strong departures from LTE (e.g., Eriksson & 
Toft 1979, Kiselman 1993, Shchukina et al. 2005). Accurate 
oxygen model atoms are available in the literature (e.g., 
Allende Prieto et al. 2003a, Shchukina et al. 2005), which allow 
to confidently perform non-LTE calculations. The importance 
of inelastic collisions with neutral H, often neglected in the 
non-LTE calculations, have been explored by Allende Prieto 
et al. (2004b) and shown to be necessary to accurately repro- 
duce the center-to-limb variation of the triplet line-profiles in 
the solar spectrum. Studies of the O i triplet are also important 
because the triplet lines are excellent probes of the physics of 
line formation (e.g., Reetz 1999). 



Using the exact same atomic model and non-LTE calcula- 
tions as those used in Allende Prieto et al. (2004b), 11 we com- 
puted non-LTE line-profiles and compared them to those ob- 
served in the spectra of BD +17 4708 (Fig.Qj}- LTE profiles 
were also computed for completeness. Non-LTE calculations 
were made both with and without including H collisions. In 
the former case, the simple approximation formula by Drawin 
(1968), enhanced by an empirical factor S h, as suggested by 
Steenbock & Holweger (1984), was adopted. 

The effect of including H collisions is evident from the 
fits in Fig. ^2 Our non-LTE calculations without H collisions, 
when forced to fit the reddest line of the triplet by tuning the 
oxygen abundance, overestimate the two strongest lines of the 
triplet. As S h is increased, the three lines are better fitted, si- 
multaneously, while reducing the size of the non-LTE correc- 
tion to the oxygen abundance at the same time. A very good 
fit to the BD +17 4708 triplet profile is found with S H = 10. 
Note, however, that in the solar case the S h = 1 model fits ac- 
curately the center-to-limb observations (Allende Prieto et al. 
2004b). Although Allende Prieto et al. did not test non-LTE 
calculations with S h = 10, such case would lead to an oxygen 
abundance in disagreement with the other oxygen abundance 
indicators. 

The triplet lines originate from the 3s 5 S° to 3p 5 P transition. 
The upper level consists of three states with very similar ener- 
gies but noticeably different transition probabilities. This im- 
plies that the depths of formation of the three lines are different. 
The bluest, strongest line is formed in an upper layer while the 
reddest, weakest line of the triplet is formed in a deeper layer. 
Since the density of neutral hydrogen decreases with depth in 
the stellar atmosphere, the collisional rates due to neutral H are 
more important in the upper layers thus thermalizing more ef- 
ficiently the 3p 5 P levels, which results in a weakening of the 
line strengths. Considering the different formation depths of 
the three lines and the fact that the 3p 5 P level becomes pro- 
gressively more thermalized with height due to H-collisions, 
the bluest line of the triplet gets more weakened than the red- 
dest line. This reasoning is consistent with what is shown in 

Fig. CD 

Adopting A = 8.70 for the Sun, we find an LTE [O/Fe] ra- 
tio of 0.87 for BD +17 4708. The best non-LTE fit to the data, 
that for S h = 10, reduce this ratio to 0.81. Oxygen abundances 
in metal-poor stars have been determined by several groups us- 
ing different lines and types of analyses without general agree- 
ment. The so-called 'oxygen abundance problem' is complex 
and still open (see, e.g., Nissen et al. 2002, Fulbright & Johnson 
2003, Melendez et al. 2006). Discrepancies regarding whether 
the [O/Fe] ratios remain constant at about +0.5 for the most 
metal-poor stars or if they increase as lower [Fe/H] values are 
reached still exist although most authors favor constant [O/Fe] 
ratios. In view of our results for BD +17 4708, we conclude 
that H-collisions are an important ingredient in the non-LTE 
computations of the triplet, not to be ignored in metal-poor 
stars. However, it is still unclear if the approximation adopted 



11 Allende Prieto et al. used both ID and 3D model atmospheres in 
their study. We remind the reader that our work is restricted to ID 
Kurucz models. 
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in this work is accurate, given that two different values of S h 
are needed to fit the spectra of the Sun and a metal-poor star. 
In any case, when non-LTE oxygen abundances are inferred 
from the triplet, the three observed line profiles must be accu- 
rately reproduced, simultaneously, by the models. At present, 
our analysis does not allow us to give a reliable estimate of the 
oxygen abundance of BD +17 4708. 



5. Conclusions 

The high accuracy with which the spectral energy distribution 
of the SDSS standard BD +17 4708 has been measured has al- 
lowed us to provide a reliable estimate of its effective tempera- 
ture. We have then used spectral line analysis to infer consistent 
parameters for this star. 

Once the degeneracy between r e ff and E(B - V) in the 
model fits to the observed flux distribution is broken by in- 
dependent estimates of E(B - V), which included a detailed 
modeling of interstellar absorption features in the observed 
spectrum, we obtain the following parameters, given here 
along with reasonable estimates of the (lcr) error bars: T e ff = 
6141 + 50 K, [Fe/H] = -1.74 ± 0.09, logg = 3.87 ± 0.08, 
E(B - V) = 0.010 + 0.003. The spectral energy distribution 
also allowed us to obtain reliable values for the bolometric flux, 
/boi = 4.89 + 0.10 x 10~ 9 erg cirT 2 s , and angular diameter, 
6 = 0.1016 + 0.0023 mas, of the star. We thus provide accurate 
(in an absolute sense) parameters for the spectrophotometric 
standard of the Sloan Digital Sky Survey. 

Compared to previous spectroscopic studies, our T e g is 
higher by about 190 K, which has a severe impact on classical 
abundance analyses of moderately metal-poor stars. For exam- 
ple, this increase ruins the good agreement between the iron 
abundance derived from Fe i and Fe n lines found when a lower 
T e g is used. Despite this, all other features on the spectrum 
(e.g., the Balmer lines or the strong Mg lines used to constrain 
the logg value) seem to be more consistent with the models 
when a high T^g is adopted. In particular, the excitation bal- 
ance of Fe i lines is satisfied with the high T e ff but not with the 
lower value. 

We also determine the mean abundance of ar-elements 
([a/Fe] = 0.37 ± 0.06). The non-LTE modeling of the permit- 
ted oxygen triplet lines should include the effect of collisions 
with neutral H to reasonably reproduce the observations but a 
better physical treatment of H collisions is needed. 
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Appendix A: Comparison with the literature 

Here we compare the atmospheric parameters derived in this 
paper for BD +17 4708 with those given in the literature 
(Tabled, with the exception of Peterson (1981), which is not a 
CCD-based paper. For the absolute iron abundance (Ap e ) com- 
parison, corrections due to the different solar iron abundances 
adopted in each study and effective temperature difference ef- 
fects (Fe i lines only) are taken into account. For the latter we 
adopt a correction of 0.2 dex per 300 K (6.67 x 10" 4 dex K" 1 , 
according to Fig. |9} when necessary. In each case, 'our' Ap e 
value refers to the abundance we would derive with the param- 
eters adopted by each author or group. 

- Rebolo et al. (1988) adopted a temperature scale similar 
to that used by Peterson (1981), which is coupled to older 
Kurucz models by photometric calibrations based essen- 
tially on synthetic photometry, but with different color cal- 
ibrations. They derive T e ff = 5890 K. Their \ogg = 4.0 was 
obtained from photometric calibrations based on Stromgren 
photometry. Only one Fei line was used to derive their 
Ap e = 5.80. We confirmed the accuracy of the E.W. given in 
their paper and for this line only we derive the exact same 
abundance when using the Rebolo et al. parameters. 

- Magain ( 1989) used a photometric calibration based on the 
IRFM (Magain 1987) to derive T eff = 5960 K. This rela- 
tively low IRFM r e ff is due to the assumption of zero red- 
dening. The logg = 3.40 derived in this work is too low 
and it does not satisfy ionization equilibrium with our line 
data when the low T e ff is adopted. Some, but not all, the gf 
values used in this work are based on a solar analysis. They 
derive Ap e = 5.75, which is about 0.05 dex lower than our 
result for the Fe i lines. 

- Axer et al. (1994) inferred their T e $ = 6100 K from model 
fitting of the Balmer lines. Note that their high logg = 4.4, 
obtained by forcing ionization equilibrium of Fe lines, is 
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reasonably expected for the temperature adopted. They ob- 
tain A Fe = 6.09 while for their T e ff our Fe i lines suggest 
Ap e = 5.90. Axer et al. noticed that due to systematic dif- 
ferences in their E.W. measurements compared to at least 
two previous studies, their abundances are probably over- 
estimated by about 0.15 dex. The S/N in this study is sig- 
nificantly lower than in most others. If we correct for this 
likely error in the E.W. measurements, then A Fe = 5.94, 
in good agreement with our result. However, the gf val- 
ues used by these authors were determined using the solar 
spectrum so the good agreement may be fortuitous. 
Spite et al. (1994) do not give enough details to make a fair 
comparison. Their T e ff = 5950 K and logg = 3.30 values 
are from the excitation and ionization balance conditions, 
but no details are given about the line list and atomic data. 
No solar A Fe is given either. 

Thevenin & Idiart (1999) used the Thevenin (1998) cat- 
alog as the source for their LTE parameters but no de- 
tails on their determination are given, except that it is a 
re-analysis of literature data. They derive non-LTE correc- 
tions to the Fe abundances, which amount to about 0.2 dex 
in the case of BD +17 4708. According to the authors, 
only the Fei lines suffer from significant deviations from 
LTE. With their non-LTE corrections, the iron abundance 
increases from A Fe = 5.71 to 5.92. Our LTE abundance for 
their r e ff is Ap e = 5.81, from both Fei and Fe n lines (within 
0.02 dex). The 0.1 dex difference in the LTE abundance is 
likely due to the use of solar gf values in their work. As it 
is discussed in 33.61 the increase in the Fe i abundance due 
to the non-LTE effects predicted by these authors worsens 
the ionization balance problem we find. 
Boesgaard et al. ( 1999) used two T e ff scales, those by King 
(1993, K93) and Carney (1983a,b; C83), which result in 
6091 K and 5956 K, respectively. The former is based on 
the modeling of the Ha line while the second is essen- 
tially the Peterson (1981) T e ff scale. The low T e ff obtained 
in this way may be due to several factors, for example: 
1) missing opacity in older Kurucz models which leads to 
overestimated UV and visible fluxes and thus lower T e g to 
match the observations, 2) metallicity effects not properly 
accounted for in the T e ff vs. color calibrations, which gen- 
erally result in low temperatures for metal-poor warm stars 
if the calibration is constructed mainly with solar metallic- 
ity stars, 3) the zero point correction to the T e ff scale. Their 
[Fe/H] values have been taken from the literature but put 
onto the same scale by using the same solar iron abundance. 
They find A Fe = 5.78 with the K93 scale and A Fe = 5.70 
with the C83 scale. These values are both 0.10 dex larger 
than ours. 

Fulbright (2000) performed a classical spectroscopic abun- 
dance analysis, deriving from the excitation equilib- 
rium of Fe i lines condition and then setting log g from the 
ionization balance condition using Fen lines. The process 
is iterative but it starts with the r e jf estimate. Given that 
log g has a smaller effect on the Fe I abundances, the result- 
ing temperatures and metallicities are mostly affected by er- 
rors in the Fe i line modeling. The complete line list used in 
this study is published and a throughout comparison can be 



made. In fact, we reproduced this analysis using the same 
models (Kurucz overshoot) and atomic data. Our resulting 
Ap e vs. E.P. relations are shown in Fig. IA.ll The gf values 
given by Fulbright have been empirically corrected so it is 
not surprising to find a smaller scatter compared to our re- 
sults (we remind the reader that we avoided this type of cor- 
rections given that they may artificially reduce the impact 
of model uncertainties). As shown in Fig. IA.ll with the T e ff 
adopted by Fulbright (6025 K), a small trend with E.P. still 
remains but it disappears for T e ff = 6190 K. In Fig. IA.ll we 
also show the slope of the A Fe vs. E.P. trends (e) as a func- 
tion of r e ff (for comparison purposes we also show our e 
values). This T e ff increase degrades the almost perfect ion- 
ization balance obtained with T e ff = 6025 K (where the Fe i 
and Fe n abundances agree within 0.03 dex compared to a 
difference of about 0.08 dex with the high T e ff). Thus, it 
seems that the small deviation from excitation balance was 
sacrificed to almost perfectly satisfy the ionization balance 
condition. Clearly, the line-to-line scatter in this analysis 
is still too large as to use the excitation/ionization balance 
conditions as good T^l \ogg indicators, i.e., there will be 
always room to alter the T e ff and log g values by changing 
the criteria for excitation and ionization balance, which are 
not going to be satisfied simultaneously due to model limi- 
tations. Note also that the E.P. coverage of the Fulbright line 
list is about 2 eV shorter than ours (at least for the lines that 
can be reasonably well analyzed). The Fulbright analysis 
suggests A Fe a 5.90, which is larger by about 0.1 dex than 
our A Fe , most likely due to the use of Kurucz overshoot in- 
stead of no-overshoot models. Despite the empirical correc- 
tions, the gf scale of Fulbright is in good agreement with 
the one used in this work. The use of the Barklem et al. 
(2000) damping constants instead of the modified Unsold 
approximation values used by Fulbright does not change 
the abundances by more than 0.01 dex given that the lines 
are relatively weak so this does not explain the excitation 
balance discrepancy. Neither does the use of no-overshoot 
models. 

- Mishenina et al. (2000) obtained - 5943 K from the 
Alonso et al. (1996) IRFM temperature scale and 6000 K 
from their fitting of models to the Ha line. Both values 
need to be corrected upwards due to reddening and prob- 
ably observational errors. They derive A Fe = 5.90, in rea- 
sonable agreement with our A Fe = 5.84 (after correcting 
for T e g and solar abundances). Their logg = 4.0 value 
was inferred by forcing ionization equilibrium of Fe lines. 
However, their mean Fe i and Fe n abundances have uncer- 
tainties of 0.11 and 0.17 dex respectively. Although they 
use the correct logg, the error bar in their derived value 
should be very large and the claimed agreement between 
Fe i and Fe n abundances is questionable. 

- Ryan et al. (2001) determined T e ff - 5983 K by using sev- 
eral color calibrations, including that by Magain (1987), 
in principle all giving essentially the same T e g. It is inter- 
esting, then, to find that even adopting E(B - V) = 0.01, 
they find a low T e g with a T e s scale consistent with that by 
Magain (1987). This is likely because Ryan et al. explicitly 
avoided the high temperatures predicted by the IRFM (see 



Ramirez et al.: The SDSS standard BD +17 4708 



17 



6.5 



•> 6.0 



5.5 



6.5 



v 6.0 



5.5 



<B o 



T „ = 6025 K 



T e „ = 6190 K 



0.0 1.0 2.0 3.0 4.0 
E.P. (eV) 



5.0 





0.04 


> 






0.02 


X 









-3 


0.00 


w 






-0.02 



II II II M_ 


J II 1 | 1 II 


M I 1 | 1 1 1 1 


1 II 1 II 


II J I 1 J I 1 J 


MM Ill 


Mini 





5900 6000 6100 6200 6300 
T e „ (K) 

Fig.A.l. Top panel: abundance of iron as a function of exci- 
tation potential for two different effective temperatures as in- 
ferred using the Fulbright (2000) line list. Open circles: Fei 
lines, filled circles: Fe n lines. The solid lines are linear fits to 
the Fe i data only. Bottom panel: slope of the Ap e vs. E.P. trend 
as a function of 7^ using the Fulbright (2000) line list (squares 
and solid line) and our line list (triangles and dotted line). One- 
sigma (lcr) errors for the slopes are also shown. 



§5.2 in Ryan et al. 1999), claiming an unphysical nature of 
the color-calibrations by Alonso et al. (1996), which have, 
however, been recently confirmed by Ramirez & Melendez 
(2005b, their §4.4). The Ryan et al. [Fe/H] and log g values 
were taken from the literature. 

- Simmerer et al. (2004) determined T e ff from the Alonso 
et al. (1996) calibrations but with E(B - V) = 0, hence 
the low T eS =i 5941 K. They inferred logg = 3.98 both 
from the Hipparcos parallax of the star and the ionization 
balance condition of iron lines, which is in good agree- 
ment with our results for the low T e «. Their iron abun- 
dance is Ap e = 5.92 while ours for the same parame- 
ters is Ap e = 5.82. The 0.1 dex difference between these 
abundances is probably due to the use of very weak lines 
(E.W. < 10 mA in most cases), which suffer from errors 
due to noise and continuum placement. Note that the gf 
values used in their work are on the same scale as ours (i.e., 
they were obtained from the same sources), but the lines are 
not the same. 

- Nissen et al. (2004) obtained T e « = 5943 K, also from the 
Alonso et al. (1996) calibrations with zero reddening (as 



inferred from the Schuster & Nissen 1989 calibration) and 
logg = 3.97 from the Hipparcos parallax. They find Ap e = 
5.89, in reasonable agreement with our Ap e — 5.82. Note 
that they used MARCS instead of Kurucz models. Most of 
their iron lines are in the blue (A < 4600 A) so it is not 
possible to make a line-by-line comparison. 

- Melendez & Ramirez (2004) used literature values for 
[Fe/H] and \ogg. Their T e ff = 6154 K is from the Ramirez 
& Melendez (2005b) IRFM temperature vs. color calibra- 
tions, using E(B - V) =* 0.02. Note that for E(B - V) = 0, 
the IRFM temperature vs. color calibrations suggest a lower 
T eff 6050 K (Fig. ED while application of the IRFM for 
this star results in 5955 K (Ramirez & Melendez 2005a). 
At E(B - V) = 0.01, however, the r e jf from the color cali- 
brations is in good agreement with that found in our work. 

- Asplund et al. (2005) determined = 6183 K from fit- 
ting of the wings of the Ho- line using essentially the same 
atomic data that we used in 33 .41 but MARCS models in- 
stead of Kurucz. The effect of using a different model at- 
mosphere is negligible given that we find a very similar Ha 
temperature ( 33.4> . Their logg value is inferred from the 
Hipparcos parallax of the star. They derive [Fe/H] = -1.51 
(Ap e = 5.99) from both Fei and Fen lines, which gave 
almost the same value for this star (but they find a mean 
systematic difference of 0.08 dex with the Fei lines giv- 
ing lower abundances in their complete sample). The use 
of the MARCS model adopted by Asplund et al. increases 
the Fe n abundance by 0.06 dex with respect to the abun- 
dance we derived with the Kurucz model but the Fe i abun- 
dance remains unchanged. Asplund et al. used only weak 
Fe i lines with gf values from O'Brian et al. (1991) and Fe n 
lines for which gf values are from Biemont et al. (1994). 
When we restrict our line list to lines with these character- 
istics, the Fe i abundance reduces by 0.03 dex while the Fe n 
abundance increases by 0.02 dex. Thus, the mean (Fei- 
Fe n) difference reduces from 0. 15 dex when using our line 
list, adopted atomic data, and a Kurucz model to 0.03 dex 
if a MARCS model and a line selection similar to that of 
Asplund et al. is made. 

In summary, the effective temperature differences between 
the one obtained in this paper and those given in the litera- 
ture can be explained by ignored reddening, limitations of older 
models, errors in the observations and basic data in the Ha line 
modeling, and the nature of the classical spectroscopic anal- 
ysis. In most cases we were able to reproduce the literature 
abundances with corrections due to solar Ap e and r e (j differ- 
ences. In the other cases, the remaining differences (0.1 dex or 
less) can be reasonably explained by the atomic data adopted 
and/or different model atmospheres. This means that the codes 
for LTE analyses are free of major bugs. Our Ap e values, with 
conservative error bars, are very robust given that our line list 
is accurate in an absolute sense 



